Effect of transport-induced charge inhomogeneity on point-contact Andreev reflection 
spectra at ferromagnet-superconductor interfaces 
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We investigate the transport properties of a ferromagnet-superconductor interface within the 
framework of a modified three-dimensional Blonder- Tinkham-Klapwijk formalism. In particular, we 
propose that charge inhomogeneity forms via two unique transport mechanisms, namely, evanescent 
Andreev reflection and evanescent quasiparticle transmission. Furthermore, we take into account the 
influence of charge inhomogeneity on the interfacial barrier potential and calculate the conductance 
as a function of bias voltage. Point-contact Andreev reflection (PCAR) spectra often show dip 
structures, large zero-bias conductance enhancement, and additional zero-bias conductance peak. 
Our results indicate that transport-induced charge inhomogeneity could be a source of all these 
anomalous characteristics of the PCAR spectra. 
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I. INTRODUCTION 



Measurement of spin polarization is an important sub- 
ject in spintronics. The point-contact Andreev reflec- 
tion (PCAR) technique has been widely used as a sim- 
ple and powerful method in the determination of con- 
duction electron spin polarization in a broad range of 
ferromagnets[ll-Q. The principle behind the spin polar- 
ization measurement by the PCAR spectroscopy is based 
on the fact that the Andreev reflection [4] probability at 
a ferromagnet-superconductor interface is limited by the 
minority spin carrier density at the Fermi level in the fer- 
romagnet. In the limit of a clean ballistic ferromagnet- 
superconductor contact, the spin polarization P is simply 
determined by the ratio G(V — 0)/Gn of the conduc- 
tances G(V = 0) and Gn at zero and high bias voltage, 
respectively, which was obtained by decomposing the cur- 
rent into a fully polarized part and a fully unpolarized 
part.[l],[2| However, for real ferromagnet-superconductor 
contacts, which are generally not in the clean limit, ac- 
curate determination of P is nontrivial and usually re- 
quires a careful analysis of the complete conductance- 
voltage (G — V) curve. [3[ Blonder et al. [5j developed 
a theory (known as the BTK model) that takes into ac- 
count the Andreev reflection, for the experimental G — V 
curves of normal metal-superconductor contacts with dif- 
ferent behaviors ranging from metallic to tunnel junction. 
Later, Strijkers et al. p] extended the BTK model to a 
ferromagnet-superconductor contact by including the ef- 
fects of the spin polarization in the ferromagnet. The 
extended BTK model of Strijkers et al. has been ap- 
plied to analyze the G — V curves of many ferromagnet- 
superconductor contacts and also to determine the spin 
polarization in the ferromagnets in the contacts. 0, [f| 
In addition to the original one-dimensional (1-D) BTK 
model .V and the extended 1-D BTK model|6], several 
extended three-dimensional (3-D) BTK models for both 
ballistic and diffusive regimes have also been proposed. 0- 



Interestingly, anomalous phenomena such as the dip 
structures 12. ITll-ll9| . the zero-bias conductance peak 



(ZBCP) ljj|20j-[22j, and the large zero-bias enhancement 
QjJ , may appear in the PCAR spectra of a wide variety 
of point contact systems. For example, Fig. 1(a) shows 
that unexpected dips emerge when the bias voltage V ap- 
proaches the superconducting energy gap A. In addition, 
a narrow peak occurs at V = (known as the ZBCP), as 
shown in Fig. 1(b) and (c). The large zero-bias enhance- 
ment may result in the normalized zero-bias conductance 
being larger than 2.0 [see Fig. 1(a)]. These character- 
istics cannot be explained by the above-mentioned ex- 
tended BTK models 0, [EI-fToj and thus other effects must 
be taken into account. 

The physical origins of the dips, ZBCP, and large zero- 
bias enhancement have been studied before. For ex- 
ample, Strijkers et al.® investigated the dip structures 
by introducing an additional proximity-induced super- 
conducting energy gap. On the other hand, intergrain 
Josephson effect [23| and Maxwell resistance [24j were also 
regarded as possible sources of both dips and large zero- 
bias enhancement. Schmidt et a/. [251] proposed formation 
of superconductor-insulator-superconductor (SIS) junc- 
tions that would result in the ZBCP. Indeed, the ge- 
ometry of a point contact may lead to the formation 
of a SIS junction. However, the ZBCP also appears 
in planar superconductor-insulator-normal metal (SIN) 
junctions [i~7|. In ref. [26j . the ZBCP was explained 
by multiple phase-coherent reflections and proximity- 
induced Josephson effect. However, this explanation is 
inapplicable for the material with large exchange split- 
ting. Furthermore, the authors did not provide any 
model to fit their experimental data. Interestingly, the 
Andreev bound states (27} have also been considered as 
a source of the ZPCP. A recent theoretical work further 
showed that an Andreev bound state could be viewed as a 
topological edge state [28|. However, an Andreev bound 
state stems from the unconventional pairing symmetry 
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FIG. 1. (Color online) The anomalous PCAR spectra from 
our own measurements at T — 4.2K (solid lines) (a-c). (a) Dip 
structures for a Nio.8iFeo.19/Pb point contact, (b) The zero- 
bias conductance peak for a Nio.8iFeo.19/Pb point contact, 
(c) The zero-bias conductance peak for a Gao.94Mno.06As/Pb 
point contact. In (d), the previously measured PCAR spec- 
trum with a broadened superconducting energy gap for a 
Fe/Pb point contact |6|, together with the fit to the extended 
BTK model[6j (the dashed line), is also shown. 



and thus would not occur in s-wave superconductors. 
Therefore, the origins of the dip structures, zero-bias en- 
hancement, and ZBCP are still not fully understood. 

In addition to the above-mentioned anomalous fea- 
tures, for some experimental data, the value of the fit- 
ted model parameters for the extended BTK model of 
Strijkers et o/.[B], e.g., A, is much larger than the well- 
known value. For example, in Fig. 1(d), the energy gap 
of Pb is 1.35 meV, whereas the fitted value is 1.93 meV. 
Such PCAR spectra which exhibit a broadened energy 
gap, are called broadened spectra. However, two plau- 
sible mechanisms: thermal effects [29l - [3^ | and a spread- 
ing resistance [33U35T ] . have been proposed to explain the 
broadened spectra. 

In this paper, we investigate the transport mechanism 
of a ferromagnet-superconductor interface and develop a 
modified BTK model to better describe the physics at 
the interface and account for anomalous PCAR spectra. 
We solve a 3-D Bogolubov de-Gennes equations, assum- 
ing no scattering along the transverse direction. We find 
that evanescent waves form for the large-angle incidence. 
We point out that these evanescent waves might lead 
to charge inhomogeneity and an electric dipole layer oc- 
curs at the interface. Furthermore, we take the effect 
of charge inhomogeneity into account by modifying the 
barrier strength Z. Interestingly, by doing that, we es- 
tablish a new 3-D BTK model which can describe the 
anomalous conductance spectra with dip structures, a 
zero-bias conductance peak, and large zero-bias enhance- 
ment. In other words, we show that all anomalous con- 
ductance spectra have the same origin: transport induced 
charge inhomogeneity. We also discuss the temperature- 
dependence of the dip structures. Futhermore, we fit the 
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FIG. 2. (Color online) (a) For the spin-up incident electron, 
the incident region could be divided into three parts. Here, 9 
is the incident angle. In region I, all wave functions are trav- 
eling waves. In region II, the wave function of an Andreev re- 
flected hole would decay. In region III, all the wave functions 
of Andreev reflected hole, hole-like, and electron-like states 
are evanescent. The boundaries of the three regions depend 
on the spin polarization (P). (b) A schematic plot of charge 
distribution near the interface due to the evanescent Andreev 
reflection and evanescent quasiparticle transmission, (c) The 
relation between contact area and incident angle. The white 
region on the left represents the ferromagnet while the gray 
(blue) region on the right indicates the superconducting tip. 
The black circle denotes the incident electron. 



experimental results using this new model including the 
effect of spreading resistance. We find that not only the 
ZBCP spectra but also unsual spectra could be fitted well 
by our model. 



II. THEORETICAL MODEL 

Wc begin with the 3-D Bogolubov de-Gennes Hamil- 
tonian. 
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Here, exchange potential h(r) = E ex Q(—x) and pair po- 
tential A(r) = Ae^O(cc), the cc-axis is normal to the 
interface. The ferromagnet and superconductor are on 
the left and right side of the interface, respectively (see 
Fig. 2). a denotes up (down) spin and p^ (pi) is 1 (-1). 
Like other BTK formalisms, we model the interface scat- 
tering by a delta- function potential U5(x). In addition, 
we define the spin polarization 



(N t v Ft - iV^UFj.) 
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where A^(A^) and vf^(vf±) are density of states and 
Fermi velocity for spin-up (down) electrons, respec- 
tively. We obtain P — E ex /Ep- Furthermore, we as- 
sume the E and \JE 2 — A 2 terms are negligible and 
hence obtain the wave vectors for electrons, holes, and 
quasiparticles in superconductor \k ea \ s» kp^/1 + p a P, 
\kha\ ~ — Po-P, \k s \ ~ hp, respectively. Here, 

k F = y/2m(E F + E ex )/h. 



A. Incident angle and Andreev reflection 

We consider allpossible incident angles which vary 
from to 7r/2@, Q and assume no scattering parallel 
to the interface (i.e. the transverse component of wave 
vector k is the same in the ferromagnet and superconduc- 
tor) 0. We then derive the eigenstates of the Bogolubov 
de-Gennes Hamiltonian. 

By matching the transverse components of wave vec- 
tors at the boundary, we find that with for the spin- 
up incident electrons, the forms of wave functions could 
be separated into three distinct regions, as shown in 
Fig. 2(a). When the angle of incidence is less 
than cos _1 [ v / 2P/(l + P)], wave functions do not de- 
cay. Conversely, when the incident angle is larger 
than cos -1 \\J2P / (1 + P)], the Andreev reflected hole be- 
comes evanescent. Increasing the incident angle above 
cos _1 [-^/P/(l + P)], the transmitted states in the super- 
conductor (i.e. electron-like and hole-like states) start to 
decay and cannot propagate. Therefore, we use the two 
boundary conditions: *&f(x = Q — ) = ^s(x = + ) and 
<F£ s I 'dx\ x=0 - - d^ F /dx\ x=0 + = 2Zk F y>(0) to match the 
wave functions for the three regions, separately. Here, 
Z is the barrier strength. Following the original BTK 
model [a], Z is defined by mU/h 2 kp. 

In region I (i.e. cos 2 9 > 2P/(1 + P)), the wave func- 
tions in the ferromagnet and superconductor read as fol- 
lows. 
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(3) 



In region II (i.e. P/(l + P) < cos 2 (6>) < 2P/(1 + P)), 
unlike region I, we take into account the evanescent hole 
state and hence ^ f is written as 



^F = 



(4) 



However, is the same as in region I. 

In region III (i.e. cos 2 9 < P/(l + P)), in addition to 
the evanescent Andreev reflected hole, the transmitted 
waves decay exponentially and hence we have to rewrite 
*s as, 



*« = c 



u 



"0 



(5) 



where uo and vo have the same definition as in ref. 

The wave vectors in the x-direction are k ex = 
k F (l + P) 1 / 2 cos6», k hx = k F [(l + P)cos 2 9 - 2P] 1 / 2 , 
and k sx = k F [(l + P) cos 2 9 ~ P] 1 / 2 . The decay pa- 
rameters ph x = kp[2P — (1 + P) cos 2 6*] 1 / 2 , and p sx — 
k F [P — (1 + P) cos 2 9} 1 / 2 . kF denotes the Fermi wave 
vector. 

As to the spin-down incident electron, the wave func- 
tions never decay regardless of how large the incident 
angle is. Therefore, the forms of wave functions are the 
same as that in region I except the wave vectors are 
different. The wave vectors in the x-direction for the 
spin-down incident electron are k ex — kp[L — P] 1 / 2 cos#, 
k hx = k F [(l - P) cos 2 9 + 2P] 1 / 2 , and k sx = k F [(l - 
P) cos 2 + P] 1/2 . 

The wave vectors of the incident electron and the An- 
dreev reflected hole are not equal. This indicates that the 
reflected hole does not retrace the path of the incident 
electron, as in the usual Andreev reflection process. As 
pointed out in ref. (36|, for the 3-D model, in the presence 
of the exchange interaction, the retro-reflectivity of An- 
dreev reflection would be broken. Recently, it was shown 
that specular Andreev reflection may hap pen in graphene 
metal-superconductor structures [371 138| . As to the de- 
caying waves at larger incident an gle, they stand for the 
states localized near the interface [39j]. After matching 
wave functions at the interface, we can obtain the am- 
plitudes of the Andreev reflection a, normal reflection b, 
electron-like transmission c, and hole-like transmission d 
and hence the probabilities for each region. 

In addition, it should be noted that evanescent 
Andreev reflection and evanescent quasiparticles tran- 
simission would also happen for a normal metal- 
superconductor interface. The evanescent Andreev 
reflection and quasiparticle transmission occur when 
fc e sin# > kh and fc e sin# > k s . When k e = kh (k e = k s ), 
k e sin 9 < kh [k e sin 9 < k s ] for any incident angle. There- 
fore, the inequality of wave vector is a necessary condition 
for the formation of the evanescent states. For higher spin 
polarization (ferromagnet), the inequality of wavectors is 
mainly caused by the exchange splitting. Hence, it is 
reasonable to neglect the terms of E and \J E 2 — A 2 for 
mathematical simplication. However, when the spin po- 
larization is close to zero, E and \JE 2 — A 2 are not neg- 
ligible compared to the exchange splitting. For P w 0, 
if we keep E and \J E 2 — A 2 in wave vectors, the wavec- 
tors of the incident electron and Andreev reflected hole 
are given by (Hk e ) 2 /2m = (hkp) 2 /2m + E, (hkh) 2 /2m = 



{hk F ) 2 /2m - E, (Hk+) 2 /2m = {Hkp) 2 /2m + VE 2 - A 2 , 
and (hk-) 2 /2m = (hk F ) 2 /2m - \J E 2 - A 2 . Here, fc+ 
and fcr are the wave vectors of electron-like and hole-like 



states, respectively. In the presence of E and y/E 2 — A 2 , 
therefore, the evanescent Andreev reflection and quasi- 
particle transmission would appear in the normal metal- 
superconductor point contacts as well. Our mathemat- 
ical formulae presented in this paper are applicable for 
the ferromagnet-superconductor interfaces only. In the 
future, we will study the evanescent transport at normal 
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metal-superconductor interfaces. 

B. Modification of barrier strength 

In the original BTK theory |5j, it was assumed that the 
applied bias voltage has no effect on the barrier poten- 
tial or Z. However, this assumption was not based on 
any physical argument and may be invalid. From the 
above calculation, we find that evanescent waves form 
for the large-angle incidence. Therefore, the likely influ- 
ence of evanescent states on the barrier potential of the 
interface should be considered. The evanescent wave rep- 
resents the spatially inhomogeneous distribution of parti- 
cles. Hole and quasiparticles carry charges. This implies 
that charge may not spread uniformly and could be local- 
ized near the interface instead. This charge inhomogene- 
ity would affect the transport of other incident electrons 
via electrostatic interaction. 

In the evanescent Andreev reflection process, the inci- 
dent electron near the interface would take another elec- 
tron away from the ferromagnet to the superconductor, 
and an evanescent hole is therefore left in the ferromag- 
net. The movement of a hole in one direction implies the 
movement of an electron in the opposite direction. The 
evanescent Andreev reflected hole does not carry any cur- 
rent. That is, the removed electron does not travel either 
after it is transferred into the superconductor. Therefore, 
via the evanescent Andreev reflection, positive charges 
would prefer to reside on the ferromagnetic side of inter- 
face. Conversely, negative charges would be localized on 
the superconducting side of the interface. Hence, an elec- 
tric dipole layer forms, as shown in Fig. 2(b). Similarly, 
electric dipole layer could also occur via the evanescent 
quasiparticle transmission because the charges of inci- 
dent electrons are transferred into the superconductor 
via the formation of decaying quasiparticles and local- 
ized positive charges appear in the ferromagnet near the 
interface. The related study on dipole layer at metal- 
superconductor interface was conducted in ref. (40| . In 
addition, as shown in Fig. 2(a), evanescent transport 
happens at larger incident angles. Furthermore, in Fig. 
2(c), an electron at the black circle may be incident at 
the ferromagnet-superconductor interface along various 
directions. However, the maximum possible incident an- 
gle is limited by contact area. For example, the contact 
area A\ is larger than A2 and hence the maximum pos- 
sible incident angle Q\ > 62- The maximum possible 
incident angle increases with contact area. Therefore, we 
can infer that the effect of electric dipole layer would be 
more striking when the contact area is larger. 

The electrostatic potential due to an electric dipole 
layer may be written in the form(4]| 

*(x)^/^ (x >.V'(^)^. (6) 

Here, -D(x') is the surface dipole moment density and 
n is the direction of dipole moment of the layer. Obvi- 



ously, the electrostatic potential decreases with the dis- 
tance between the observer and dipole layer, and becomes 
quite strong when the observer is very close to the dipole 
layer. The electric dipole layer due to charge transfer at 
the ferromagnet-superconductor interface should display 
the similar behavior and hence has a maximum in the 
vicinity of the interface. Furthermore, taking into ac- 
count screening effect, the range of the potential could 
be finite. In metals, the screening length is about a 
few Angstroms. Therefore, it would be reasonable to 
model the screened short-range electrostatic potential as 
a 5-function. Using the definition of the charge den- 
sity Q = e(|/| 2 - \g\ 2 ) in ref. @, we find that the 
charge density of an evanescent hole is proportional to 
a* a. Therefore, the barrier potential due to the contri- 
bution of an evanescent Andreev reflection is linear in 
a* a. However, the strength of barrier potential is dif- 
ficult to determine, because the (5-function is unphysical 
in reality [29| and merely for the purpose of mathematical 
simplification. Therefore, we set the strength of barrier 
potential equal to 1/2(0-* a), where U2 is an undetermined 
constant. We use the similar method to model the barrier 
potential due to evanescent quasiparticles transmission. 
Furthermore, using Eq. (6), we can find that the electro- 
static potential energy for incident electrons due to these 
evanescent states is negative. Because a in the region 
III is found to be 0, the barrier potential is modified as 
U = Ui - U 2 (a*a) n - U 3 {c*c - d*d) ln (u^u - v$v ). 
We can derive the modified barrier strength 

Z = Z 1 -Z2{a*a) II -Z 3 (c*c-d*d) III (u* ) u -VQV ). (7) 

Here, Z\ is related to the scattering due to the contact po- 
tential, impurity, and tunneling barrier at the interface, 
as the previous BTK defined. Z2 is the barrier strength 
that the evanescent Andreev reflection contributes. Z 3 
is caused by the Coulomb potential due to the evanes- 
cent quasiparticles transmission. We assume that the 
localized charges affect only the transport of the incident 
spin-up electrons in region I and the incident spin-down 
electrons. In other words, we neglect the influence the 
localized states in regions II and III on the transport of 
each other. The barrier potential for the spin-up incident 
electrons in regions II and III is not altered by the charge 
inhomogeneity. Therefore, we use this modified Z for the 
incident spin-up electrons in region I and the incident 
spin-down electrons. For the spin- up incident electrons 
in regions II and III, the barrier strength is equal to Z\. 
The [a*a)ii and (c*c — d*d)in(uQU Q — VqVq) are shown 
in Table. I. 



C. Calculation of conductance 

Finally, by averaging over incident angles between 
and 7r/2, as in ref. @, we obtain the average Andreev 
reflection Aj-(A^) and normal reflection B^(B^) probabil- 
ities for the spin-up (down) electron. Using these proba- 
bility coefficients, we can calculate the conductance of a 
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TABLE I. The analytical forms of (a*a)ii, where a is the complex coefficient in Eq. (4), and (c*c — d* d) 1 1 1 (u^uo — vqVq) in 
Eq. (5). Here, ug = 1 - = \{l + y/1 - A 2 /E 2 ) and A = P/(l + P) cos 2 9. 
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point contact as in the BTK model 0. The conductance 
is given by 



G(V) = N^Avpt 
+ N ie 2 Av Fi 



df(E-eV) 

dE 
df(E-eV) 

dE 



(1 + Af - B t )dE 

(1 + A± - BJdE 
(8) 



Here. 



both 



Nf/Ni and Vp\/vp^ are equal to 
Furthermore, we normalize the con- 



P)/(l-P) 

ductance by Gat, which is equal to G(V)| e v»A- 

In order to explain the spectra broadening effect, we 
also take into account the extra resistance factor T, 
which results from the spreading resistance. In the pres- 
ence of spreading resistance, the normalized conductance 
Gfs/Gfn is modified as 



g fs \ = i + r 
Gfn I {Gfs/Gfn)^ 1 +r 

with r = R s / RpN, 



(9) 



where Rfn = ^/Gfn represents the resistance at large 
bias voltage, while R s denotes the spreading resistance. 
We use this formalism to fit the experimental data and 
the results are presented in Sec. IV. 



III. EXPERIMENTS 

We fabricated some permalloy (Nio.8iFeo.19) thin films 
on Si substrates by DC magnetron sputter system. The 
thickness of samples are about 100 nm. The bulk resis- 
tivity is ^287 (nfim). To perform the PCAR measure- 
ments, a differential screw was used to bring a fine tip 
onto the sample surface. Fine Pb tips fabricated by me- 
chanical grinding were pressed onto the sample repeat- 
edly to ensure good contact between the tips and the 
samples. To minimize thermal fluctuation, conductance- 
voltage curves were taken by the differential technique^ 
in a helium cryostat for 4.2 K measurements. The con- 
ductance curves were normalized by the high bias values. 
Each sample was measured 50 times on different contact 
condition. 



The contact size was estimated from the quasiclassical 
formula of contact resistance, applicable to both ballistic 
and diffusive transport (26l l43j 



R= (1 



3ird 2 1 \ d 2d 



(10) 



where d is the diameter of the contact, p is the resistivity 
of the material, ^(im/d) is a factor in the order of unity. 
From this contact resistance formula we can calculate 
the contact diameter in the range of about 2 ~ 19 nm. 
Selected PCAR spectra with dip structures and ZBCP 
from our permalloy thin films are displayed in Fig. 1(a) 
and (b), respectively. Some of the previously measured 
dilute magnetic semiconductor samples |34j (Fig.l (c)), 
which have much larger resistivities, can also be analyzed 
with this model. 



IV. RESULTS AND DISCUSSION 

A. Effects of barrier strengths 

In Fig. 3, the calculated conductance spectra under 
three different barrier strengths are displayed in order 
to examine the influence of the barrier strength. We 
find from Fig. 3(a) that when the effect of charge in- 
homogeneity are not taken into account (i.e., Z2 = 
and Z 3 — 0), the normalized conductance at low bias 
decreases with increasing barrier strength Z\ and two 
peaks appear symmetrically about zero bias, like the pre- 
vious BTK models @, Q. This was already attributed 
to the suppression of Andreev reflection by the barrier 
strengthQ. In Fig. 3(b), it is shown how the evanes- 
cent Andreev reflection [through barrier strength (i.e. 
Zi 7^ 0)] affect the conductance spectra. Clearly, the 
low-bias conductance would increase initially, with in- 
creasing Zii as shown by the curves with Zi = 0.1 and 
Z-2 = 0.4, but, however, would decrease upon further in- 
creasing Z2. From Eq. (7), we can infer that Zi would re- 
duce the effective barrier strength Z and hence increases 
the low-bias conductance. However, if Z2 is large, the 
effective barrier strength becomes negative. The nega- 
tive potential barrier implies a potential well rather than 
a potential barrier. When the effective barrier strength 
becomes negative, the depth of the well increases with in- 
creasing Zi and hence the conductance would be reduced 
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FIG. 3. (Color online) (a) The conductance spectra for 
different origin barrier strength at fixed spin polarization 
P = 0.3 when Z2 and Z3 are not applied, (b) The effect of 
evanescent-hole-induced barrier strength Z2 on conductance 
spectra, (c) The effect of evanescent-quasiparticle-induced 
barrier strength Z3 on conductance spectra. The parameters 
P, A, and T in (b) and (c) are the same as (a). 



by Z 2 . Moreover, when the conductance increases with 
Z 2 (i.e. the curves for Z 2 = 0, 0.1, and 0.4), the dip at the 
zero-bias would diminish gradually and eventually a sin- 
gle peak occurs at zero bias (see the curve for Z 2 — 0.4). 
The shape of conductance spectra for Z 2 = 0.4 is similar 
to that of the curve in the clean regime (i.e. the curve for 
Z\ = 0, Z 2 = 0, Z3 = 0, shown in Fig. 3(a)). However, in 
comparison with the curves with no barrier, the height 
of the peak for Z 2 = 0.4 is higher. The heights of peaks 
for Z 2 = 0.4 and for no barrier arc 1.81 and 1.67, respec- 
tively. In addition, we find that the peak for Z 2 = 0.4 
is sharper. We also notice that Z 2 can generate double 
dips, like the curves for Z 2 = 0.8,0.9, 1.0 and 1.1. 

In Fig. 3(c), the effect of the evanescent-quasiparticle- 
induced barrier strength Z3 is shown. In contrast with 
Z 2 , Z3 would reduce the low-bias conductance first and 
then raise it. It is clear from Table I that the coefficient 
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FIG. 4. The simulated anomalous conductance spectra with 
(a) dip structures, (b), (c) zero-bias conductance peak, and 
(d) both dip structures and large zero-bias enhancement. 



(c*c - d*d)iu(uQUo - VqVq) ^ for \E\ > A and the 
(a*a)u approaches zero at high energy. Therefore, G^v 
is insensitive to Z 2 but may be influenced by Z3. Similar 
to Z 2 , for smaller Z3, Z3 reduces the height of the barrier 
at high bias and hence GV increases and the low-bias 
normalized conductance decreases. For sufficiently large 
Z3 , a barrier is transformed into a well and Gjq decreases 
with Z3. Hence, the normalized conductance could be 
enhanced by increasing Z3. Furthermore, Z3 could raise 
zero-bias conductance above 2 (see the curve for Z3 = 
2.6). 

It is worth mentioning that, under the effect of charge 
inhomogeneity, the zero-bias normalized conductance 
does not reach its maximum in the clean limit. Indeed, 
the transmission of electrons into the superconductor 
should be obstructed by a barrier (or a well) especially 
upon Andreev reflection, as shown in Fig. 5 of ref. 0. 
Furthermore, at high bias, electrons propagate into the 
superconductor mainly through electron- and hole- like 
transmission and hence Gat is independent of Andreev 
reflection probability. Therefore, the zero-bias conduc- 
tance G(0)/Gn must decrease if Z is a constant. When 
Z 2 and Z3 are included, the total barrier strength varies 
with voltage. Electrons may feel a higher barrier or a 
deeper well at high bias and therefore normalized con- 
ductance could increase in the presence of a barrier or a 
well. 



B. Anomalous spectra 

From the preceding section, we may conclude that the 
evanescent- hole-induced barrier strength Z 2 could lead to 
zero-bias enhancement, sharp central peak, occurrence 
of dips while the evanescent-quasiparticle-induced bar- 
rier strength Z3 could substantially enhance the zero-bias 
conductance. By applying Z 2 and Z 3l we can create the 
anomalous PCAR spectra for different spin polarizations, 
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FIG. 5. (Color online) The temperature-dependence of dip structures, (a) The calculated conductance spectra for different 
temperatures T at fixed superconducting energy gap A. (b) The calculated conductance spectra for different superconducting 
energy gaps A at fixed temperature T. (c) The measured point-contact Andreev reflection spectra for a Fe/Nb point contact 
at different temperatures from ref. [24|. All curves in (a)-(c) except the top curve are shifted downward for clarity. The 
parameters P, Z\, Z%, and Zz in (b) are the same as (a). The numbers marked in (a), (b), and (c) are the positions of dips. 



which are very similar to experimental observations, as 
shown in Fig. 4. In other words, all the observed anoma- 
lous characteristics can be reproduced in our model with 
suitable values of the various barrier strength parame- 
ters. 

Fig. 4(a) shows the anomalous spectra with the dip 
structures. In experiments, the dips are usually located 
at eV > A@, H3. The dips occur at V = ±1.69 mV in 
the present case. In the presence of the spreading resis- 
tance, the dips would move farther away from the central 
peak. This is consistent with the experimental observa- 
tions. Fig. 4(b) and (c) show two different types of the 
ZBCP spectra. In Fig. 4(b), a higher conductance peak 
is located at zero bias and two small peaks appear sym- 
metrically on both side of it, similar to the experimental 
results in ref. [26]. In Fig. 4(c), a zero-bias shallow 
hump appears at tunneling-like Andreev reflection spec- 
tra, similar to the experimental results in refs. (l2l. |20|- 
[22l ] . The two ZBCP characteristics can both be obtained 
by adjusting Z2, as shown in the curves for Z2 = 0.8 
and Z 2 — 0.9 of Fig. 3(b). Our results indicate that 
both dip structures and a ZBCP are due to the evanes- 
cent Andreev reflection. In ref. [12], it was showed that 
the dip structures in the Andreev spectra [like Fig. 4(a)] 
could be transformed into the ZBCP [like Fig. 4(b)] by 
varying contact area, which is consistent with our results. 
Adjusting contact area could alter the properties of in- 
terface, i.e., Z%, Z2, and Z3 and hence different types of 
the PCAR spectra would occur. 

Fig. 4(d) shows both the dip structures and large 
zero-bias enhancement. In experiments, the dip struc- 
tures are often accompanied by the large zero-bias 
enhancement 0, [23] . Our model can reproduce both the 
dip structures and large zero-bias enhancement. In addi- 
tion, we find from Fig. 3(b) that the dips appear when Z2 



reduces the zero-bias conductance. This indicates that 
the unusual zero-bias enhancement results from Z3. In 
experiments, the dip structures, a ZBCP, and large zero- 
bias enhancement usually happen at large contact area 
(i.e. low resistance). Charge inhomgeneity is also pro- 
nounced for large contact area, as seen in Fig. 3(c). This 
supports the validity of our model. However, it should 
be noted that the increase of conductance with contact 
area is due to the increase of electron-transport chan- 
nels instead of Z2 or Z3. In the preceding subsection, it 
was illustrated how Z 3 reduces Gm and G(0) under the 
assumption of constant electron-transport channels. In 
principle, the charge inhomogeneity could also happen 
for the planar geometry. Therefore, our model could also 
be extended to the planar junctions. 

C. Temperature dependence of dip structures 

Let us now examine how the dip structures vary with 
temperature. The superconducting energy gap decreases 
with increasing temperaturefiij . We plot the normalized 
conductance for different values of temperature and en- 
ergy gap in Fig. 5(a) and Fig. 5(b), respectively. Fig. 
5(a) indicates that the dips would be reduced gradually 
and also move away from the central peak as T increases. 
Fig. 5(b) shows that, with decreasing A, the dips would 
diminish gradually, too, but move towards the central 
peak. We can infer that the dips would disappear even- 
tually when the temperature is raised. Moreover, from 
Fig. 2-2 in ref. [44|, near T = 0, the energy gap is 
almost insensitive to T. If T is larger, A would de- 
crease with T more rapidly. When T is increased up 
to the critical temperature T c , A drops to zero. In other 
words, we predict that, with the temperature rising, the 
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FIG. 6. (Color online) The conductance spectra for 
Nio.8iFeo.ig/Pb point contacts measured at 4.2K (black cir- 
cles) with contact resistance of (a) R = 3.80 SI and (b) 
R = 8.83 Q,. The solid (red) and dashed (blue) lines are 
the fits using our model (see Sec. II) and the extended BTK 
model[(|, respectively. The fit value of temperature is fixed 
at 4.2 K. 



D. Fitting experimental spectra 

We performed the least squared fit to selected experi- 
mental data using our model and compared with the ex- 
tended BTK model. @ In addition to the P, Z\, Z 2 , and 
Z$, which are necessary in our model, A, T, and the extra 
resistance factor T in Eq. (9) can also be used as fitting 
parameters to the experimental spectra. To minimize the 
fitting parameters, we chose to fit the measurements ob- 
tained from samples immersed in liquid Helium and fixed 
the temperature T — 4.2 K. At first, superconducting en- 
ergy gap A for Pb was fixed at 1.23 meV and the T was 
fixed at zero. If the fitting result was not satisfactory, we 
would turn on A and/or F as fitting parameters. 

We find that, for small contact resistance (R = 3.80 
f2), a zero-bias conductance peak could appear in the ex- 
perimental curve [Fig. 6(a)]. However, for large contact 
resistance (i? = 8.83 f2), we do not observe the ZBCP and 
an ordinary PCAR spectrum occurs, as shown in Fig. 
6(b). We also find from Fig. 6(a) that, in contrast to 
the extended BTK model d] which could not reproduce 
the shoulder structures, our model could fit the ZBCP 
spectrum very well. Fig. 6(b) shows that our model also 
fits usual spectra well. Although the usual spectra could 
be fitted by the extended BTK modcl[6| well, the re- 
sultant parameter P = 0.24 is significantly smaller than 
that from the previous measurements pi I45rj47j . The spin 
polarization of permalloy has been determined by both 
PCARdlHl and spin-polarized tunneling [46l. l47j exper- 
iments, and was found to be in the range of 0.3-0.5. In 
our model, the fitted value of P ~ 0.45 agrees well with 
these previous measurements. This indicates that, in our 
measurements, charge inhomogeneity may happen and, 
therefore, the extended BTK model could not explain 
the experimental data in Fig. 6(b). 



position of a dip should increase if temperature is low 
enough. Conversely, at higher temperatures, due to the 
rapid decrease of A, the position of a dip could decrease 
with increasing temperature. Fig. 5(c) shows the mea- 
sured temperature dependence of the dips for a Fe/Nb 
point contact [13], which should correspond to the the- 
oretical results in Fig. 5(b). Many other experimen- 
tal works have shown the same temperature dependence 
of dip structures for a wide variety of combinations be- 
tween metals and superconductors [ill, H3- Fur- 
thermore, it was pointed out in the preceding subsection 
that dip structures and the ZBCP have the same physi- 
cal origin. Therefore, indeed, the two local minima in the 
ZBCP spectra would be equivalent to the dip structures. 
Kastalsky et al. [13] observed ZBCP spectra in the tem- 
perature range of T = 0.5 — 1.7K. Their data showed the 
two local minima in the conductance spectra moved to 
higher bias as temperature increased. Therefore, we infer 
that the dips may also move away from the central peak, 
as shown in Fig. 5(a). 



V. CONCLUSIONS 

In conclusion, we have developed a modified 3-D BTK 
model by taking into account the effect of transport in- 
duced charge inhomogeneity on the barrier potential at 
ferromagnet-superconductor interfaces. To this end, we 
add two new parameters, namely, the evanescent hole- 
induced barrier strength Zi and evanescent quasiparticlc- 
induced barrier strength Z3 to describe the modified bar- 
rier strength. We find that Z2 would give rise to the 
dip structures and zero-bias conductance peak, while, on 
the other hand, Z3 would enhance the zero-bias conduc- 
tance by decreasing the normal-state conductance Gjv- 
We have, therefore, offered a possible explanation for the 
anomalous conductance spectra which have been long 
standing problems in the field of PCAR spectroscopy. 
Furthermore, the results of our model fits to the exper- 
imental data have a good agreement with the measured 
ZBCP spectra as well as normal conductance spectra. 
Nevertheless, the standard <5-function form adopted here 
for the barrier potential at point contacts could be an 
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oversimplication. Further theoretical and experimental 
works are certainly needed to fully understand these im- 
portant issues in the PCAR spectroscopy. 
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